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Abstract 

Conjectured links between the distribution of values taken by the characteristic polynomials 
of random orthogonal matrices and that for certain families of L-functions at the centre of the 
critical strip are used to motivate a series of conjectures concerning the value-distribution of 
the Fourier coefficients of half-integral weight modular forms related to these L-functions. Our 
conjectures may be viewed as being analogous to the Sato- Tate conjecture for integral weight 
modular forms. Numerical evidence is presented in support of them. 



1 Introduction 



The limiting value distribution of the Fourier coefficients of integral weight modular forms is given 
by the celebrated Sato- Tate conjecture. Our purpose here is to identify the implications of some 
recent conjectures concerning the value-distribution of certain families of L-functions at the centre 
of the critical strip for the distribution of the Fourier coefficients of related half-integral weight 
modular forms. 

These conjectures concern the relationship between properties of L-functions and random 
matrix theory. It was conjectured by Montgomery that correlations between the zeros of the 
Riemann zeta function on the scale of the mean zero separation coincide with those between the 
phases of the eigenvalues of unitary matrices, chosen at random, uniformly with respect to Haar 
measure on the unitary group, in the limit of large matrix size. This is supported both by theoretical 
[221 1271 121 El and extensive numerical [231 125j evidence. Katz and Sarnak ^^1 then generalized 
the connection by suggesting that statistical properties of the zeros within various families of L- 
functions, computed by averaging over a given family, coincide with those of the eigenvalues of 
random matrices from the various classical compact groups, the particular group being determined 
by the family in question. Based on these ideas, a link was proposed in [T7] between the leading 
order asymptotics of the value distribution of the Riemann zeta function on its critical line and 
that of the characteristic polynomials of random unitary matrices, giving, for example, an explicit 
conjecture for the leading order asymptotics of the moments of the zeta function. This approach 
was then extended to relate the value distribution of L-functions, in families, at the centre of the 
critical strip, to that of the characteristic polynomials of matrices from the various classical compact 
groups 101 El . It has also recently been extended to include all lower order terms in the asymptotics 
0. (For other related results, see jl4l I15j .^ These developments have recently been reviewed in 

Our strategy here is to combine these random-matrix-theory-inspired conjectures for the 
value distribution of L-functions with formulae due to Waldspurger |H21, Shimura Kohnen 
and Zagier 1201, and Baruch and Mao jlj which relate the values taken by L-functions associated 
with elliptic curves at the centre of the critical strip to the Fourier coefficients of certain half- 
integral weight modular forms. This approach was first outlined in [HI, where its implications 
for the vanishing of L-functions were examined. Here we use the same ideas in order to develop 
conjectures for the value-distribution of the Fourier coefficients. 

The outline of this paper is as follows. In Section |21 we give a brief overview of the rela- 
tionship between modular forms and L-functions associated with elliptic curves. In Section we 
review the results from Random Matrix Theory that we shall need, and some conjectures they 
suggest for the value-distribution of L-functions associated with elliptic curves. In Section we 
combine results from Section|21and SectionOto motivate conjectures for the moments of the Fourier 
coefficients, the value distribution of the logarithm of the Fourier coefficients, and the distribution 
of the coefficients themselves. In Section jSl we review the implications, ffi'st outlined in ,9^, for 
the vanishing of L-functions at the centre of the critical strip. Finally, in Section El we describe 
numerical experiments whose results support our conjectures. 
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2 Modular forms and L-functions 



We review in this section the connection between the quadratic twists of modular L-functions and 
the Fourier coefficients of half- integral weight modular forms. This is central to the motivation 
underlying the conjectures we make in subsequent sections. In order to be concrete, we concentrate 
on the case of quadratic twists of L-functions associated with elliptic curves. 

Let Le{s) be the L- function associated with an elliptic curve E over Q with Dirichlet series 
and Euler product given by 

oo 

n=l p|A pfA 

= n^p(l//), K(s)>3/2, (2) 

p 

with A the discriminant of E, and Op = p + 1 — ^E{¥p), where ^E{¥p) denotes the number of 
points on E when regarded over Fp. It is a consequence of the recently solved Taniyama-Shimura 
conjecture |33j j31j [Ij that Le{s) has analytic continuation to C and satisfies a functional equation 
of the form 

27r , , ^ , , / 27r ^ ''"^ 



T{s)LEis)=WE[^] r{2 -s)LEi2-s), (3) 



/QJ ^ ' ' - \^ 
where Q is the conductor of the elliptic curve E and we = il- 

We let Xd{fi) = (^) for a fundamental discriminant d, where (^) is the Kronecker symbol. 
The twisted L-function 

LE{s,Xd) = }_^ (4) 

n=l 

is the L-function of the elliptic curve Ed-, the quadratic twist of E by d. If (d, Q) = 1, then Le{s, Xd) 
satisfies the functional equation 

27r / 27r \ ''"^ 

r{s)LE{s,Xd)=Xd{-Q)wEi^=-j:) r{2-s)LE{2-s,Xd). (5) 



We shall be interested in the case when WEXd{—Q) = 1) since otherwise LE{l,Xd) is trivially equal 
to zero. 

We now come to a key point. The L-functions above are related to half- integer weight 
modular forms via formulae due to Waldspurger ,32_, Shimura j3Uj . Kohnen and Zagier and 
Baruch and Mao 1 . One must distinguish between positive and negative d, and one must also 
sort them according to various residue classes. This has been worked out explicitly for thousands 
of examples by Rodriguez- Villegas and Tornaria in the case that Q is squarefree and they kindly 
supplied a database of such forms to the authors 24; . 

Specifically, for Q squarefree, assume that d < and that Xdip) = — «p for every p \ Q 
(op = ±1 for p I (5 when Q is squarefree). Notice that such d are restricted to p\Q iip ~ 

podd 

residue classes mod Q or AQ, depending on whether Q is odd or even. Then, for such d, 

LE{l,Xd) = KE j==-, (6) 
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where C£;(|d|) S TL are the Fourier coefficients of a weight 3/2 modular form, and where ke depends 
on E. 

For d > 0, the work of Baruch, Mao, Rodriguez-Villegas, and Tornaria ^ gives the 
relevant weight 3/2 form for Q prime and d's satisfying Xd(Q) = ^Q- Oris has the same relation 
as above, but with a different proportionality constant. Some examples are listed in Sectional In 
either case, given a coefficient of the weight 3/2 form, the constant ke can be evaluated either by 
comparison with the Birch and Swinnerton-Dyer conjecture, or by numerically computing Le{1, Xd) 
as a series involving the exponential function. 

Our strategy will be to write down conjectures for the value distribution of the Fourier 
coefficients c{\d\) by coupling the connection (jB)) with conjectures motivated by random matrix 
theory for the value distribution of LE{l,Xd)- 

3 L-functions and random matrices 

It was conjectured by Montgomery 122' that the zeros of the Riemann zeta function are distributed 
statistically like the eigenvalues of random hermitian (self-adjoint) matrices, or, equivalently, like 
the phases of the eigenvalues of random unitary matrices. This extends to the zeros of any given 
principal L-function |2Zj. It was conjectured by Katz and Sarnak that the distribution of zeros 
defined by averaging over families of L-functions with the height up the critical line fixed coincides 
with the distribution of the phases of the eigenvalues of matrices from one of the classical compact 
groups, depending on the family in question. 

These ideas motivated the conjecture jl7j that, asymptotically, the moments of the Riemann 
zeta function (or any other principal L-function) averaged high on its critical line coincide, up to a 
simple arithmetical factor, with the moments of the characteristic polynomials of random unitary 
matrices. This suggestion was then extended to relate the moments of families of L-functions at the 
centre of the critical strip to those of the characteristic polynomials of matrices from the various 
classical compact groups jSllTHj. 

For any elliptic curve E it is conjectured that the family of even- functional equation 
quadratic twists 

= {Le{s, Xd) ■■ WEXd{-Q) = +1} (7) 

is orthogonal. Specifically, this family conjecturally has symmetry type . Thus the value dis- 
tribution of LE{l,Xd) should be related to that of the characteristic polynomials of matrices in 
SO{2N), at the spectral symmetry point, with N ~ log(|d|). 

For an orthogonal matrix A, the characteristic polynomial may be defined by 

Z {A, 9) = det (^I - Ae-'^y (8) 

The eigenvalues of A form complex conjugate pairs e^*^", and so the symmetry point is at = 0. 
The moments of Z{A, 0) are defined by averaging over A with respect to normalized Haar measure 
for S0{2N), dA: 

Mo{N,s)= [ |det(/-^) I'd^. (9) 

JsO{2N) 
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It was shown in JH] that for Res > — 1/2 
and that as — > oo 

Mo(iV,s)~<?.(0+)iV^(^-i)/2, (11) 

with 

^ ' VG(i + 2s)r(i + s) 

where G is Barnes' G-function: 

oo 

G{z + 1) = (2^)^/2 exp(-((7 + 1)^2 + ^)/2) TT (l + exp(-z + zV2n). (13) 



n=l 



It follows from the fact that G(l) = 1 and G{z + 1) = T{z)G{z) that 

fc— 1 

5.(0+) = 2'=('^+i)/2n A, (14) 
i=i 

for integer k. Note that the right-hand side of ()1U|) has a meromorphic continuation to the whole 
complex s-plane. 

For integer k, Mo{N, k) is a polynomial of order k{k — l)/2: 

Mo(A^,fc)= |2'=('=+i)/^n-|-) n (^ + ^)- (15) 



o<j<i<fc-i 



It can also be written in terms of a multiple contour integral i?}, a form that will be useful later 
for comparison with L-functions: 

X n (17) 

l<e<m<k llj=l^j 

Here A(xi, . . . ,x„) = ni<i<j<n(^j ~ ^^^^ contours enclose the poles at zero. 

The value distributions of the characteristic polynomial and its logarithm can be written 
down directly using (|1()|) . Let P/v(t) denote the probability density function associated with the 
value distribution of |Z(A, 0)|, i.e. 

meas.{^ G S0{2N) : a < |det(/ - A)\ < P} = / PN{t)dt. (18) 

Ja 

Then 

f c+ioo 



PN{t) = ^J Mo{N,s)t-'ds, (19) 
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for any c > 0. 

The asymptotics of -P/v(t) as t ^ comes from the pole of Mo{N,s) at s = —1/2. The 
residue there is 

^ ' T(N)^^T{j-l 2)T{i + N -2, 2) ^ ' 
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which, as ^ oo, is given asymptoticahy by 



(Ar)llr(j-l/2)rO- + iV-3/2) 
icahy by 

h(N) ~ 2~^I^G(\l2y-^I^N^I^. (21) 



Thus 

P/v(t) ~ h{N)t-^l'^ (22) 

as t — > 0. 

Importantly for us, one may deduce a central limit theorem for log|Z(yl, 0)| from equa- 
tion (UHJ (see HH]): 



A^^oo VlogTV V^Ja 2 



The above results, proved for the characteristic polynomials of random matrices, suggest 
the following conjectures for L£;(l,Xd) (these are special cases of those made in jTS], but with the 
number theoretical details worked out explicitly). 

Recall that we are assuming that Q is squarefree. Let {d,Q) = 1. For LE(s,Xd) to have an 
even functional equation, one needs WEXd{—Q) = 1- This imposes a condition, in the case that Q 
is odd, on d mod Q, and, in the case that Q is even on d mod AQ (4(5 because Xd(2) is periodic 
with period 8). Let 

Q \Q if Q is odd and squarefree ^^^^ 
1 4(5 if Q is even and squarefree. 

Next, we focus our attention on a subset of the d's according to certain residue classes mod Q. We 
let 

S-{X) = S^{X) = {-X <d<0; Xd{p) = -dp for ah p \ Q] (25) 

i.e. the set of negative fundamental discriminants d up to —X, but restricted according to a 
condition on d mod Q. Let 

and 

M^(^,s) = ^z^ E LEil^XdT. (27) 

' ^ des-{x) 

i.e. the sth moment of LE{l,Xd)- 

For elliptic curves E whose conductor Q is prime, we also look at the set of positive funda- 
mental discriminants 

S+{X) = S+{X)={0<d< X- XdiQ) = ag} (28) 
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and define W^{X,t), M^{X, s) as in the negative case. 
Note that 

W^iX,t) = ^ y^_^^ M^{X,s)t-'ds. (29) 

The central conjecture is that as X ^ oo 

M^{X,s) ~ A^{s)Mo{logX,s) (30) 

where 

A^{s) = n (1 - P-Y'-''^' (^) [I + \ W/pY + £,(-l/p)^)) (31) 
P\Q 

is an arithmetical factor that depends on E. A heuristic for A^{s) is given in [S]. The following 
conjectures are then motivated by this. 

First we consider the moments of elliptic curve L-functions. 
Conjecture 1 For M^{X,s) defined as in i27\ ), 

Second, for large X, as t ^ we have 

WiiX,t) A^i-^)hilogX)t-'/^ (33) 
and thus by scaling t so that t = 0((log X)"'^), for 7 > 1, we are led to 

Conjecture 2 // logX/(logX) —^OasX—^ 00, then 



hm y3H4^Wl{XJ{logX)y) = By-^' (34) 
(logX)8 

where 

5 = 2-7/8G(l/2)7r-i/^^±(-l/2) (35) 

In addition, we have 
Conjecture 3 



lim 



Vlog log|d| 



j exp(--)dt. 



/27r 

We take log -L£;(l, Xd) + ^loglog|d| to lie outside the interval if LE{l,Xd) = 



(36) 
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The above conjecture is similar to central limit theorems for the Riemann zeta function and 
other L-functions usually attributed to Selberg [29} 12811 . An analogous conjecture is made in JH] 
for quadratic Dirichlet L-functions. 

Further, we have a conjecture, closely related to the preceding one, for the distribution of 
the full range of values of LE{l,Xd)- 

Conjecture 4 



for fixed < a < p. 

Finally, it is discussed in [8; how (|16|) leads to conjectures for mean values of L-functions, 
not just at leading order as in 1)32(1 . but including all terms in the expansion down to the constant 
term. 

Conjecture 5 With k an integer 



1 



{d e S^{X);a< (7b^L£;(l,Xd))Vi°gi°sM < P} 



(37) 





(38) 



as X — > oo, where is a polynomial of degree k{k — l)/2 given by the k-fold residue 




■k, (39) 



where the contours above enclose the poles at zero and 



FH^u...,z,) = Ai{z,,...,z,)]j(^!^^(^Yy n C(l + ^. + ^,) (40) 





with, for p\Q, 



-{'^iy'{r'2 (n^- (tW) {t^))) ■ (^^> 



and, for p \ Q, 




(43) 
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4 Conjectures relating to the value distribution of the Fourier 
coefficients of half-integral weight forms 

Our goal now is to use the conjectures listed at the end of the previous section for the value 
distribution of the L-functions associated with elliptic curves to motivate conjectures for the value 
distribution of the Fourier coefficients of half- integral weight forms. These follow straightforwardly 
from the connection 

In each case, we let c{\d\) refer to the coefficients of the corresponding weight 3/2 modular 
form, as in ©, and refer to the corresponding proportionality constant. 

Our first conjecture then follows from (|.S6j) : 



Conjecture 6 (central limit conjecture 

1 



lim 



X^oo \S^{X)\ 

_ / exp(--)* 



Y^log log\d\ 



(44) 



We take 2log c{\d\) — ^log{\d\) + ^loglog|(i| to lie outside the interval if c{\d\) = 



This leads directly to a conjecture for the appropriately normalized distribution of the 
coefficients themselves, which is analogous to ()c{7() . 



Conjecture 7 



1 



lim 

xZo |5±(X)| 

1 r^i _i 



\/2vr J a ^ 



{de S^{X);a < 



(45) 



for fixed < a < f3. 



Our third conjecture follows from © and 1)32(1 : 



Conjecture 8 (moment conjecture) 



(.o.,.;...-./. 154)1 ^ (^-^ 



Further, we have the conjecture following from @ and (|34() . 



8 



Conjecture 9 // logX/(logX) ^ as X ^ oo, then 



lim 



4/(logX) 1 



{deS^{Xy,^^<f{\ogX)y} 



By- 



(47) 



where 



B = 2-'^'^G{l/2)^-^I^A^{-l/2) 



(48) 



Lastly, we have the analogue of (p8|) 



Conjecture 10 With k an integer, and summing over fundamental discriminants, 



c{\d\) 



2k 



X 



T± (log(t)) (it + 0(X-2+^) 



(49) 



as X ^ oo, where Tf is the polynomial of degree k{k — l)/2 given in Conjecture\^ 



The numerical evidence that supports these conjectures is amassed in Section El 



5 Frequency of vanishing of L-functions 



Examining the frequency of L£;(l, Xd) = as d varies, as well as the order of the zeros, has particular 
significance in the context of the conjecture of Birch and Swinnerton-Dyer which says that Le{s) 
has a zero at s = with order exactly equal to the rank of the elliptic curve E. Random matrix 
theory appears to have a role in predicting the frequency of such zeros. We have argued in [U] that 
since © implies a discretisation of the values of L£;(l, Xd) and 

rwi{X,t)dt (50) 
Jo 

is the probability that Le{1, Xd) has a value of a or smaller, then the combination of ()29() and 1)30(1 
suggest the following. 



Conjecture 11 There is a constant > such that 

ties±{x) 1 

- c^X-'/\logXf/\ (51) 

I d I prime 



This conjecture first appeared in 'W , but was stated with c^ > 0. However, it became clear in 
preparing numerics for this paper that can equal zero, and an arithmetic explanation has been 
given for one of the examples in our data by Delaunay • 
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Here the fundamental discriminants have been restricted to prime values to avoid extra two 
divisibility issues placed on the coefficients c(|d|). The constant remains somewhat mysterious. 

The random matrix model suggests that it should be proportional to A^{—l/2)^J k^. However, 
when one attempts to apply the random matrix model to the problem of the discrete values taken 
on by the c{\d\) one ignores subtle arithmetic. It seems, from numerical experiments, that one needs 
to take into account further correction factors that depend on the size of the torsion subgroup of 
E, but this is still not understood. See Section 6. 

Let q \ Q he a fixed prime. Another conjecture that follows from this approach concerns 
sorting the d's for which LE{l,Xd) = according to residue classes mod q, according to whether 
Xdiq) = 1 or -1. Let 

S dGS±{X) 1 
i-i5{l>Xd)=0 

R^iX) = (52) 

2^ d6S±(X) J- 
i-B(l.Xd)=0 
Xd('?) = -1 



Conjecture 12 0/ Let 



Then, for q \ Q, 



1/2 



lim Rf(X) = R,. (54) 

A^OO 



We believe this conjecture to hold even if we allow d to range over different sets of discriminants, 
such as \d\ restricted to primes (though in the latter case we must be sure to rule out there being 
no vanishings at all due to arithmetic reasons). 

A more precise conjecture given in jJJ incorporates the next term. The lower terms do 
depend on whether we are looking at S~^{X) as opposed to S~{X). We require some notation. Let 
p be prime. For p \ Q set 

m = PHp) = i'"^^''^^^' + " + (55) 
I log(p)/(l — p) in the — case 



and for p^ Q set 



where 



flip) = 1 - ap/p + l/p 
f2{p) = 1 + ap/p+l/p. 



Next, let 



where 7 is Euler's constant. 
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Conjecture 13 ,17/ For q\Q 

+ l + aqj g + X-i{q) 



where 



<7 = -log(XQV2/(2vr))-l (59) 



These conjectures are supported by numerical evidence that will be described in the next section. 



6 Numerical Experiments 

We present numerics for 2398 elliptic curves and millions of quadratic twists for each curve 
{\d\ < 10^) confirming the aforementioned conjectures. To test these conjectures we used the 
relation To make our examples explicit, we list in Table 01 relevant data for 26 of the 2398 
elliptic curves examined. The remaining data may be obtained from the L-function database of one 
of the authors (SSJ. We used as the starting point for our computations a database of Tornaria and 
Rodriguez- Villegas |21] that lists, for thousands of elliptic curves, the relevant ternary quadratic 
forms. 

Each entry in this table contains the following data 



name 


[ai, 02, as, 04, ae] 


K 


Number of d 


number ternary forms 


relevant residue classes \d 


mo 


d Q 




linear combination 


ternary forms 



The name we use for an elliptic curve E corresponds to Cremona's table ^0], but with an 
extra subscript, either 'i' or 'r', standing for imaginary or real respectively, to specify whether we 
are looking at quadratic twists LE{s,Xd) with d < or d > 0. The naming convention used by 
Cremona includes the conductor Q in the name for the elliptic curve. So, for example, the first 
entry has name 11 Ai which is the elliptic curve of conductor 11. The 'i' indicates that we are 
looking at quadratic twists of 11^ with d < 0. 

[ai,a2,a3, 04, uq] refers to the coefficients defining the equation of E, 

+ aixy + 03^ = + a2X^ + a^x + a^. (60) 

For each curve, our data consists of values of L£;(l, Xd) with \d\ < 10^, and d restricted to p\Q {{p~ 

podd 

l)/2) residue classes mod Q or 4(5, depending on whether Q is odd or even. For imaginary twists, 
all the Tornaria-Rodriguez examples have Q squarefree, and in the case of real twists, all their 
examples have Q prime. The relevant residue classes and relevant modulus is listed in the second 
line of each entry, and the number of such d up to 10^ is given by the 'Number of d\ 

Our L{l,Xd) values are expressed in terms of the coefficients of a weight 3/2 modular 
form of level 4Q in the case of imaginary quadratic twists, d < 0, and of level 4Q^ in the case 
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of real quadratic twists, d > 0. This weight 3/2 form is expressed as a hnear combination of 
theta series attached to positive definite ternary quadratic forms. The number of forms, say r, is 
given in the last part of the first line, while the linear combination, [ai, . . . , Or], and ternary forms 
occupy the last two lines of each entry. Each ternary form is specified as a sextuple of integers 
P = [Pi, P2, ^3, Pi, (3^, Pe]. The ternary form is 

fpix, y, z) = Pix^ + P2y^ + P^z^ + P^yz + P^xz + P^xy (61) 

and the theta associated series is given by 

ep(w) = w^^^'^'y''^ (62) 

Given this data, one defines 

00 r 

J24n)w^ = Yaj9p^{w). (63) 

n=l j=l 

Then, for fundamental discriminants lying in the relevant residue classes in the table (and d < or 
d > 1 according to the name of the entry), one has 

LE{l,X<i) = ^c{\d\f/\d\'/''- (64) 

Note, our values of k and c{\d\) differ slightly from the values given in in that, in our tables, we 
ignored the value of c(l), for example in the case that d > 0, and normalized the Oj's so that the 
gcd of all the c(n)'s, for n 7^ 1, is equal to one, absorbing if necesary an extra square factor into k. 
For example, while j21j gives for the curve llAr a value of k = .25384186 . . ., we list a value of k 
that is 5^ = 25 times as big, but give c(n)'s which are 1/5 as large. 

The most comprehensive test carried out jiSj for moments involved looking at conjecture 5 
for millions of values of L^;^^^ (1, Xd), with d < and \d\ = 1,3,4,5,9 mod 11. The coefficients, 
from 8^, of the polynomials in conjecture 5 are given by Tabled In checking conjecture 5, we 
compared 

E ^iiA.(^,Xd)' (65) 

-850000<d<0 
|d| = l,3,4,5,9 mod 11 

to 

E T, (logddD). (66) 

-850000<d<0 
ld| = l,3,4,5,9 mod 11 

This comparison is depicted in Table El 

The first part of Conjecture 1 is, for integral moments, a weaker form of conjecture 5, while 
conjectures 8 and 10 follow from 1 and 5. 

In Figure n we depict the numerical value distribution of LiiA-{l,Xd), for fundamental 
discriminants —85,000,000 < d < 0, \d\ = 1,3,4,5,9 mod 11, compared to PN{t) obtained by 
taking the inverse Mellin transform of (|10|) with A'^ = 20. Because we are neglecting the arithmetic 
factor in computing the density, a slight cheat was used to get a better fit. The histograms were 
rescaled by a constant along both axis until the histogram displayed matched up nicely with the 
solid curve. Considering that we are compensating for leaving out the arithemtic factor in such a 
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r 




M2) 


/r-(3) 


/r(4) 





1.2353 


.3834 


.00804 


.0000058 


1 




1.850 


.209 


.000444 


2 






1.57 


.0132 


3 






2.85 


.1919 


4 








1.381 


5 








4.41 


6 








4.3 



Table 1: Coefficients of T~(a;) = fo{k)x^(''-^y^ + foj. k=l,2,3,4. 



k 


(El) 


m 


ratio 


1 


14628043.5 


14628305. 


0.99998 


2 


100242348.8 


100263216. 


0.9998 


3 


1584067116.8 


1587623419. 


0.998 


4 


41674900434.9 


41989559937. 


0.993 



Table 2: Moments of Lny^. versus their conjectured values, for fundamental discriminants 
-85,000,000 < d < 0, \d\ = 1,3,4,5,9 mod 11, and k = 1,...,4. The data agree with our 
conjecture to the accuracy to which we have computed the moment polynomials Tr. 



naive way, it is a bit surprising how nicely the two fit. The main point we wish to make is that the 
histogram does exhibit behaviour near the origin in support of part 2 of Conjecture 1. 

In Figure [21 we verify the central limit theorem described in Conjecture 3. The first picture 
shows the value distributions, superimposed, of 



logLj;(l,Xd) + ^loglog|d| 
Y^loglog|d| 



\d\ < 10^ (67) 



for the 26 curves described in Tables |31 The second plot depicts the average value distribution of 
the 26. These are compared against the standard Gaussian predicted in the conjecture and also 
against the density function associated to the value distribution of 



log|Z(^,0)| + llogiV 

with N = 20. This density is given by 



where 



PN{g(t))g'(t) (69) 
m = (70) 

and is shown in J18j to tend to the standard Gaussian as — > oo. To get a better fit to the 
numeric value distribution, one would also need to incorporate the arithmetic factor. In the limit, 
this factor has no effect, but the convergence to the limit is extremely slow. The variability of the 
arithmetic factor explains why the 26 densities superimposed in Figure |21 don't fall exactly on top 
of one another. To get a better fit to the average, we set p{t) = Pjq{g{t))g' [t] and plot instead 
ap{at), with a = 1.21 chosen so that the density function visually lines up with the data. 
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Figure 1: Value distribution of Le^, (1, Xd), with -85000000 < d < 0, d = 2,6,7,8, 10 mod 11. 



Conjecture 12 is verified numerically in the top plot in Figure 01 which compares, for the 
first one hundred elliptic curves E in our database, and the sets S^{X), the predicted value of Rq 
to the actual value R^{X), with X = 10^. 

The horizontal axis is q. For each q and each of the one hundred elliptic curves E we plot 
R^{X) — Rq, with X = 10^, and q < 3571. For each q on the horizontal there are 100 points 
corresponding to the 100 values, one for each elliptic curve, of R^{X) — Rq. We see the values 
fluctuating about zero, most of the time agreeing to within about .02. The convergence in X is 
predicted from secondary terms to be logarithmically slow and one gets a better fit by including 
more terms as predicted in conjecture 13. 

This is depicted in the second plot of Figure |21 which shows the difference 

R^{X) - ( 1±1Z<\ (71) 

again with X = W^, q < 3571, and the same one hundred elliptic curves E. We see an improvement 
to the first plot which uses just the main term. 

This improvement is emphasized in Figure0]which compares the distribution of R^{X) — Rq 
for all our 2398 elliptic curves, X = 10®, q < 3571, versus the distribution of (f7T|) . The latter has 
smaller variance. These distributions are not Gaussian. They depict the remainder of R^{X) 
compared to the first and second conjectured approximations. There are yet further lower terms 
and these are given by complicated sums involving the Dirichlet coefficients of Le{s). 
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-4 -2 2 4 

Figure 2: Value distribution of [log Le{1, Xd) + \ log log \d\) / y^log log \ d\ compared to the random 
matrix theory counterpart, PN{g{t))g' {t) with = 20 (rescaled as explained in the text), and its 
limit the standard Gaussian. In the first picture we superimpose the value distributions for the 26 
curves described in Table 01 The second picture shows the average value distribution of the 26. 
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Figure 3: A plot 11 for one hundered elliptic curves of R^{X) — Rq, top plot, and of (f7T|l , bottom 
plot, for 2 < g < 3571, X = 10^. 
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Figure 4: Distribution first approximation v.s. second approximation for ratio of vanishings 



To verify conjecture 11 we computed the l.h.s. of l|51|) divided by 



^±(-l/2)J4X-V4(logX)3/8 



(72) 



for our 2398 curves (545 of these have c| = and these are omitted) and X = 100000, 200000, ... 10*^. 
The resulting values are depicted in Figure IHl We also display in Figure |Slthe same data for a subset 
of 55 of these curves (those in our database with 7^ and whose conductor have leading digits 
11). The graphs displayed in these two figures appear relatively flat. 

To measure how flat these graphs are, we show in Figure [7| the distribution of the slopes of 
the graphs in Figure El measured by sampling each one at X = 5 x 10^ and X = 10^. Most have 
slopes that are of size less than 10^^'^, and this suggests that the power of X, namely 3/4, in our 
conjecture is not off by more that .01, and that the power of log(X), namely 3/8 is not off by more 
than .1. The mean does appear slightly to the right of 0, occuring at .5 x 10"^'^, but this might be 
the result of using a limited number of curves and also perhaps due to lower order terms. 

Next we sort the graphs in Figure 1^1 by their rightmost values at X = 10^, and, in Figure 
ISl plot these values against the order of the torsion subgroup of the corresponding elliptic curves. 
The plot shows that curves with trivial torsion tend to have smaller constants, followed by curves 
of torsion size equal to 7 or 5 or 3, then 2 and 1 again, etc. We do not yet have an explanation for 
this phenomenon and do not know how to incorporate it into our model. It seems j26] that to nail 
down the constant one would need to incorporate into the model Delaunay's heuristics for Tate- 
Shavarevich groups However, it appears that for primes p dividing the order of the torsion 

subgroup, the probability that c{\d\) is divisible by p deviates in a way we do not yet understand 
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from a prediction made by Delaunay for the probability that the order of the Tate-Shavarevich 
group is divisible by a given prime. Delaunay's predictions are for all elliptic curves sorted by 
conductor and here we are examining a skinny set of elliptic curves, namely quadratic twists of a 
fixed elliptic curve. 
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Figure 8: The rightmost values of the graphs in El plotted against the order of the torsion subgroups 
of the corresponding elliptic curves. We see, for example, that the smallest constants tend to go 
along with the curves whose torsion is trivial. The few curves at the far right with torsion size equal 
to 2 have c{\d\) divisible by 2 when \d\ / 4, and hence, because the corresponding discretization is 
twice as large, the c(|(i|)'s are four times as likely to vanish. 
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11 A, 


[0, -1, 1, -10, -20] 


2.91763323388 


13931691 


2 


1, 3, 4, 5, 9 mod 11 


[1/2, - 


1/2] 








[3, 15, 15, -14, -2, -2], [4, 11, 12, 0, -4, 0] 


11 Ar 


[0,-1,1,-10,-20] 


6.34604652140 


13931649 


6 


1, 3, 4, 5, 9 mod 11 


[1/10, 


-1/10, 3/10, -3/10, -2/10, 2/10] 






[1, 44, 132, -44, 0, 0], [4, 12, 121, 0, 0, -4], [4, 33, 45, -22, 


-4, 0], [5, 9, 124, -8, -4, -2], [5, 36, 36, 28, 4, 4], [16, 16, 25, 


-4, -4, 


-12] 








14A, 


[1, 0, 1, 4, -6] 


5.30196495873 


4432803 


2 


15, 23, 39 mod 56 


[1/4, - 


1/4] 








[4, 15, 15, 2, 4, 4], [7, 8, 16, 


-8, 0, 0] 






15v4, 


[1, 1, 1,-10,-10] 


3.19248444426 


4749434 


2 


2, 8 mod 15 


[1/4, - 


1/4] 








[3, 20, 20, -20, 0, 0], [8, 8, 15, 0, 0, -4] 


llA, 


[1,-1, 1,-1,-14] 


2.74573911809 


14353828 


2 


3, 5, 6, 7, 10, 11, 12, 14 mod 17 


[1/2, - 


1/2] 








[3, 23, 23, -22, -2, -2], [7, 11, 20, -8, -4, -6] 


19A, 


[0, 1, 1, -9, -15] 


4.12709239172 


14438275 


2 


1, 4, 5, 6, 7, 9, 11, 16, 17 mod 19 


[1/2, - 


1/2] 








[4, 19, 20, 0, -4, 0], [7, 11, 23,-10, -6,-2] 


19Ar 


[0, 1, 1, -9, -15] 


4.07927920046 


14438248 


12 


1, 4, 5, 6, 7, 9, 11, 16, 17 mod 19 


[-1/6, 1/6, -1/6, 1/6, -1/3, - 


1/3, -1/3, 1/3, 1/3, 1/3, 1/3, 


-1/3] 


[1, 76, 380, -76, 0, 0], [4, 20, 361, 0, 0, -4], [4, 77, 96, 40, 4, 4], [5, 16, 365, 16, 2, 4], [5, 61, 92, 16, 4, 2], [5, 76, 92, 


-76, -4, 0], [9, 44, 77, 28, 6, 


8], [16, 24, 77, 20, 8, 4], [17, 44, 44, 12, 16, 16], [20, 24, 73, 4, 8, 20], [20, 36, 45, 20, 16, 


12], [25, 36, 36, -4, -16, -16] 








21A, 


[1, 0, 0, -4, -1] 


3.82197956150 


4986931 


2 


10, 13, 19 mod 21 


[1/4, - 


1/4] 








[3, 28, 28, -28, 0, 0], [7, 12, 24, -12, 0, 0] 


26ylj 


[1, 0, 1, -5, -8] 


3.47934348343 


4704178 


2 


7, 15, 31, 47, 63, 71 mod 104 


[-1/2, 1/2] 


[7, 15, 31, -14, -6, -2], [8, 15, 28, 4, 8, 8] 




[1, -1, 1, -3, 3] 


1.80405719338 


4704185 


2 


3, 27, 35, 43, 51, 75 mod 104 


[-1/2, 1/2] 


[3, 35, 35, -34, -2, -2], [4, 27, 27, 2, 4, 4] 


SOA, 


[1, 0, 1, 1, 2] 


18.5342737810 


1583137 


2 


31, 79 mod 120 


[1/8, - 


1/8] 








[4, 31, 31, 2, 4, 4], [15, 16, 16, -8, 0, 0] 



Table 3: Some of the ternary forms from |24j . 
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33Ai I [1, 1, 0, -11, 0] I 2.74463335747 | 5224398n~3 



5, 14, 20, 23, 26 mod 33 



[1/4, 1/4, -1/2] 



[3, 44, 44, -44, 0, 0], [11, 12, 36, -12, 0, 0], [15, 20, 20, -4, -12, -12] 



MAj [1, 0, 0, -3, 1] 7.45670022989 4784604 2 



19, 35, 43, 59, 67, 83, 115, 123 mod 136 



[1/4, -1/4] 



[4, 35, 35, 2, 4, 4], [8, 19, 36, 4, 8, 8] 



35A, [0, 1, 1, 9, 1] 



2.20504427610 5540991 



1, 4, 9, 11, 16, 29 mod 35 



-1/2. 1/2] 



[4. .Sf). -M). 0. -4. ()]. [11. 15. .S9. -10. -(i. -10] 



[0, 1, 1, -23, -50] 7.07044268094 14798224 2 



2, 5, 6, 8, 13, 14, 15, 17, 18, 19, 20, 22, 23, 24, 29, 31, 32, 35 mod 37 



[1/2, -1/2] 



[8, 19, 39, 2, 8, 4], [15, 20, 23, -8, -14, -4] 



AliAr [0, 1, 1, 0, 0] 



10.9373790599 14852746 22 



2, 3, 5, 7, 8, 12, 18, 19, 20, 22, 26, 27, 28, 29, 30, 32, 33, 34, 37, 39, 42 mod 43 

[1/4, 1/4, -1/4, -1/4, 1/4, 1/4, -1/4, 1/4, -1/4, -1/4, -1/4, -1/4, -1/4, 1/4, 1/4, 1/4, -1/4, 1/4, 1/4, 1/4, -1/4, -1/4] 
[5, 69, 929, -34, -2, -2], [5, 241, 276, 104, 4, 2], [8, 65, 624, 44, 4, 4], [8, 108, 389, 88, 8, 4], [12, 29, 932, -28, -4, -4], 
[12, 73, 373, 30, 4, 8], [20, 29, 621, 22, 12, 16], [20, 61, 277, 54, 12, 8], [28, 32, 377, 28, 16, 12], [28, 33, 376, 28, 12, 
16], [29, 77, 157, 14, 18, 26], [29, 77, 161, 46, 10, 26], [32, 69, 157, 26, 12, 24], [37, 45, 237, -2, -26, -34], [37, 89, 104, 
44, 16, 10], [37, 93, 104, -60, -16, -2], [45, 77, 113, 62, 42, 10], [48, 77, 104, -56, -4, -32], [48, 77, 108, -48, -28, -32], 
[61, 80, 89, 76, 42, 16], [69, 76, 77, -4, -58, -32], [69, 77, 80, -44, -8, -58] 



67Ai [0, 1, 1, -12, -21] 6.05993680291 14974655 3 



1, 4, 6, 9, 10, 14, 15, 16, 17, 19, 21, 22, 23, 24, 25, 26, 29, 33, 35, 36, 37, 39, 40, 47, 49, 54, 55, 56, 59, 60, 62, 64, 65 

mod 67 
[1/2, -1, 1/2] 

[4, 67, 68, 0, -4, 0], [15, 36, 39, -16, -14, -4], [16, 19, 71, 6, 16, 12] 



67Ar [0, 1, 1, -12, -21] 1.27377003655 14974660 70 



1, 4, 6, 9, 10, 14, 15, 16, 17, 19, 21, 22, 23, 24, 25, 26, 29, 33, 35, 36, 37, 39, 40, 47, 49, 54, 55, 56, 59, 60, 62, 64, 65 
mod 67 

[-1/2, 1/2, 1/2, -1, -1/2, -1, -1, 1, -1, -1, -1, -1, 1, 1, -1, 1, -1, -1, 1, -1, -1, 1, 1, -1, 1, 1, -1, -1, 1, 1, -1, -1, 1, -1, -1, 1, 

1, -1, -1, 1, 1, -1, 1, 1, 1, -1, 1, 1, 1, 1, -1, 1, -1, 1, -1, 1, 1, -1, 1, 1, -1, 1, -1, -1, -1, -1, 1, 1, -1, 1] 

[1, 268, 4556, -268, 0, 0], [4, 68, 4489, 0, 0, -4], [4, 269, 1140, 136, 4, 4], [9, 269, 508, 92, 8, 6], [16, 17, 4493, 2, 16, 4], 
[17, 237, 332, -172, -8, -6], [17, 268, 332, -268, -8, 0], [21, 77, 753, -38, -2, -6], [21, 217, 308, 192, 12, 2], [21, 220, 308, 
196, 12, 16], [24, 157, 328, 48, 20, 8], [24, 236, 261, -196, -20, -12], [25, 248, 248, 228, 12, 12], [29, 56, 760, -44, -16, 
-8], [29, 89, 504, 20, 24, 26], [29, 205, 216, 88, 20, 14], [29, 205, 224, 124, 16, 14], [29, 224, 224, 180, 16, 16], [33, 144, 
292, -92, -32, -28], [33, 173, 236, 124, 8, 18], [33, 173, 237, -118, -14, -18], [36, 68, 509, 48, 20, 12], [36, 173, 216, 120, 
4, 16], [37, 77, 464, 68, 8, 26], [37, 116, 293, -64, -22, -4], [37, 188, 189, 84, 10, 32], [37, 189, 189, 110, 10, 10], [40, 
77, 449, 2, 4, 40], [40, 77, 457, 18, 36, 40], [40, 149, 237, -130, -20, -16], [49, 104, 260, -92, -24, -4], [60, 77, 277, -62, 
-4, -16], [60, 84, 277, 64, 4, 44], [60, 89, 285, -42, -44, -56], [64, 129, 160, 4, 36, 44], [65, 68, 301, -20, -2, -40], [65, 68, 
304, -12, -28, -40], [65, 140, 140, 12, 32, 32], [65, 140, 160, -84, -44, -32], [65, 149, 149, 30, 58, 58], [65, 149, 157, 110, 
34, 58], [68, 77, 237, 10, 12, 24], [68, 84, 265, 8, 16, 68], [68, 96, 217, -28, -8, -60], [68, 96, 237, 92, 12, 60], [68, 140, 
149, 72, 44, 52], [73, 77, 240, 32, 20, 46], [73, 92, 181, 12, 18, 8], [73, 132, 132, -4, -32, -32], [76, 77, 272, 48, 12, 72], 
[76, 93, 173, -6, -8, -20], [76, 93, 188, 16, 68, 20], [77, 84, 193, 40, 22, 12], [77, 89, 220, -8, -36, -70], [77, 129, 157, 70, 
74, 66], [77, 132, 148, 92, 52, 56], [77, 132, 157, 100, 74, 56], [81, 93, 196, 4, 16, 74], [84, 116, 153, 88, 32, 60], [84, 
116, 169, -48, -80, -60], [88, 92, 173, 76, 44, 36], [88, 96, 173, 80, 44, 52], [89, 132, 132, -4, -72, -72], [89, 132, 148, -92, 
-48, -72], [92, 92, 181, -12, -12, -84], [93, 93, 173, -6, -6, -82], [93, 93, 193, -46, -46, -82], [96, 96, 157, -16, -16, -76], 
[100, 104, 149, 64, 96, 44], [104, 104, 149, -64, -64, -60] 
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I [1, -1, 0, 4, -3] I 2.79278430294 



14992818 4 



5, 7, 10, 11, 13, 14, 15, 17, 20, 21, 22, 26, 28, 29, 30, 31, 33, 34, 39, 40, 42, 43, 44, 45, 47, 51, 52, 53, 56, 58, 59, 60, 

62, 63, 66, 68 mod 73 

[-1/2, 1/2, -1, 1] 

[7, 43, 84, -40, -4, -6], [11, 28, 80, 28, 4, 8], [15, 39, 40, 20, 8, 2], [20, 31, 44, -28, -4, -12] 



79Ar [1, 1, 1, -2, 0] 11.9016007052 



15008174 40 



3, 6, 7, 12, 14, 15, 17, 24, 27, 28, 29, 30, 33, 34, 35, 37, 39, 41, 43, 47, 48, 53, 54, 56, 57, 58, 59, 60, 61, 63, 66, 68, 

69, 70, 71, 74, 75, 77, 78 mod 79 

[-1/2, -1/2, 1/2, 1/2, -1/2, -1/2, -1/2, 1/2, -1/2, -1/2, -1/2, 1/2, 1/2, -1/2, 1/2, 1/2, 1/2, -1/2, 1/2, -1/2, -1/2, -1/2, 

-1/2, -1/2, 1/2, 1/2, 1/2, 1/2, -1/2, 1/2, 1/2, -1/2, 1/2, 1/2, 1/2, -1/2, -1/2, 1/2, -1/2, 1/2] 

[12, 369, 449, 54, 8, 4], [17, 133, 896, -12, -16, -14], [17, 224, 521, 20, 6, 8], [17, 264, 449, -12, -14, -16], [24, 185, 461, 
106, 4, 8], [28, 192, 385, -112, -12, -4], [29, 33, 2088, -8, -28, -6], [33, 116, 544, 92, 32, 12], [37, 153, 376, 8, 36, 34], 
[41, 56, 1044, 4, 24, 40], [41, 116, 425, -4, -38, -8], [41, 185, 265, -38, -22, -2], [41, 185, 276, -84, -32, -2], [41, 216, 
236, 20, 28, 36], [41, 217, 229, -18, -30, -14], [48, 112, 377, 28, 40, 4], [48, 193, 224, -32, -36, -20], [53, 201, 233, -182, 
-10, -30], [56, 96, 377, -44, -32, -4], [60, 157, 216, -20, -4, -32], [61, 96, 372, -4, -32, -44], [68, 113, 284, 104, 12, 20], 
[68, 172, 185, -92, -32, -4], [69, 137, 217, 10, 22, 34], [69, 161, 193, -74, -38, -14], [77, 113, 249, -62, -46, -26], [77, 137, 
201, -38, -50, -22], [85, 113, 233, 82, 46, 42], [93, 145, 161, 70, 22, 46], [96, 136, 197, -132, -60, -20], [96, 140, 185, 
-128, -16, -44], [108, 120, 161, -52, -36, -4], [108, 145, 161, 70, 36, 104], [112, 116, 185, 100, 76, 12], [113, 116, 185, 
-100, -14, -56], [113, 161, 165, 138, 2, 110], [116, 137, 185, 134, 100, 96], [132, 137, 149, 98, 64, 108], [132, 153, 156, 
-44, -124, -100], [140, 145, 153, 134, 132, 64] 



89S, [1, 1, 0, 4, 5] 



2.18489393577 



15029294 7 



3, 6, 7, 12, 13, 14, 15, 19, 23, 24, 26, 27, 28, 29, 30, 31, 33, 35, 37, 38, 41, 43, 46, 48, 51, 52, 54, 56, 58, 59, 60, 61, 

62, 63, 65, 66, 70, 74, 75, 76, 77, 82, 83, 86 mod 89 

[1/2, -1/2, 1/2, 1/2, 1/2, -1/2, -1] 

[3, 119, 119, -118, -2, -2], [7, 51, 103, -50, -6, -2], [12, 31, 92, 4, 12, 8], [15, 24, 95, 24, 2, 4], [15, 27, 96, -20, -8, -14], 

[19, 23, 95, -14, -10, -18], [23, 31, 48, 16, 12, 2] 



109Ai [1, -1, 0, -8, -7] 5.94280424076 



15060017 3 



2, 6, 8, 10, 11, 13, 14, 17, 18, 19, 23, 24, 30, 32, 33, 37, 39, 40, 41, 42, 44, 47, 50, 51, 52, 53, 54, 55, 56, 57, 58, 59, 

62, 65, 67, 68, 69, 70, 72, 76, 77, 79, 85, 86, 90, 91, 92, 95, 96, 98, 99, 101, 103, 107 mod 109 

[1/2, 1/2, -1] 

[11, 40, 119, 40, 2, 4], [19, 23, 119, -22, -18, -2], [24, 39, 56, 4, 12, 16] 



llSAj [1, 1, 1, 3, -4] 2.85781203904 15064917 7 



3, 5, 6, 10, 12, 17, 19, 20, 21, 23, 24, 27, 29, 33, 34, 35, 37, 38, 39, 40, 42, 43, 45, 46, 47, 48, 54, 55, 58, 59, 65, 66, 

67, 68, 70, 71, 73, 74, 75, 76, 78, 79, 80, 84, 86, 89, 90, 92, 93, 94, 96, 101, 103, 107, 108, 110 mod 113 

[-1/2, 1/2, -1/2, -1, -1/2, 1, 1] 

[3, 151, 151, -150, -2, -2], [12, 39, 116, 4, 12, 8], [19, 24, 119, 24, 2, 4], [20, 47, 68, -44, -4, -12], [23, 24, 119, 24, 10, 
20], [23, 40, 59, 20, 2, 8], [35, 39, 47, -10, -34, -6] 



139Ai [1, 1, 0, -3, -4] 5.80133204474 



15089693 5 



1, 4, 5, 6, 7, 9, 11, 13, 16, 20, 24, 25, 28, 29, 30, 31, 34, 35, 36, 37, 38, 41, 42, 44, 45, 46, 47, 49, 51, 52, 54, 55, 57, 
63, 64, 65, 66, 67, 69, 71, 77, 78, 79, 80, 81, 83, 86, 89, 91, 96, 99, 100, 106, 107, 112, 113, 116, 117, 118, 120, 121, 

122, 124, 125, 127, 129, 131, 136, 137 mod 139 

[-1/2, 1/2, 1/2, 1/2,-1] 

[7, 80, 159, 80, 2, 4], [11, 52, 152, 52, 4, 8], [20, 28, 139, 0, 0, -4], [20, 31, 144, 8, 20, 16], [24, 47, 71, 2, 12, 8] 
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179A, I [0, 0, 1, -1, -1] I 5.10909732904 | 1511372477 



1, 3, 4, 5, 9, 12, 13, 14, 15, 16, 17, 19, 20, 22, 25, 27, 29, 31, 36, 39, 42, 43, 45, 46, 47, 48, 49, 51, 52, 56, 57, 59, 60, 

61, 64, 65, 66, 67, 68, 70, 74, 75, 76, 77, 80, 81, 82, 83, 85, 87, 88, 89, 93, 95, 100, 101, 106, 107, 108, 110, 116, 117, 
121, 124, 125, 126, 129, 135, 138, 139, 141, 142, 144, 145, 146, 147, 149, 151, 153, 155, 156, 158, 161, 168, 169, 171, 

172, 173, 177 mod 179 

[1/2, 1/2, 1/2, 1/2,-1/2,-1/2,-1] 



[4, 179, 180, 0, -4, 0], [15, 48, 191, 48, 2, 4], [15, 51, 192, -44, 
-10], [20, 39, 184, 8, 20, 16], [39, 56, 76, -52, -20, -12] 



-14], [16, 47, 183, 6, 16, 12], [19, 39, 191, -34, -14, 



233A, [1, 0, 1, 0, 11] 1.63933561519 15133226 13 



3, 5, 6, 10, 11, 12, 17, 20, 21, 22, 24, 27, 34, 35, 39, 40, 41, 42, 43, 44, 45, 47, 48, 53, 54, 57, 59, 61, 65, 67, 68, 69, 
70, 73, 75, 77, 78, 79, 80, 82, 83, 84, 86, 87, 88, 90, 93, 94, 95, 96, 97, 99, 103, 106, 108, 111, 114, 115, 118, 119, 122, 
125, 127, 130, 134, 136, 137, 138, 139, 140, 143, 145, 146, 147, 149, 150, 151, 153, 154, 155, 156, 158, 160, 163, 164, 
165, 166, 168, 172, 174, 176, 179, 180, 185, 186, 188, 189, 190, 191, 192, 193, 194, 198, 199, 206, 209, 211, 212, 213, 

216, 221, 222, 223, 227, 228, 230 mod 233 

[1/2, -1/2, 1/2, 1, 1/2, 3/2, -3/2, 1, -1, -3, 1, 1, -1] 

[3, 311, 311, -310, -2, -2], [11, 87, 255, -82, -6, -10], [12, 79, 236, 4, 12, 8], [20, 95, 140, -92, -4, -12], [24, 39, 239, 2, 
24, 4], [24, 43, 239, 10, 24, 20], [27, 39, 244, -28, -16, -22], [35, 80, 84, 28, 24, 4], [39, 80, 96, -68, -8, -36], [40, 47, 119, 
2, 20, 8], [44, 68, 87, -36, -20, -28], [47, 68, 87, 36, 38, 40], [48, 59, 79, 2, 16, 12] 
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